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Abstract 

We examine some new DNLS-like equations that arise when considering 
strongly-coupled electron-vibration systems, where the local oscillator po- 
tential is anharmonic. In particular, we focus on a single, rather general 
nonlinear vibrational impurity and determine its bound state (s) and its dy- 
namical selftrapping properties. 
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The standard Discrete Nonlinear Schrodinger (DNLS) equation is a paradig 
matic equation that describes wave propagation in nonlinear and disordered 
systems [T]], polaron formation and dynamics in deformed media0, power- 
switching in coupled nonlinear fibers and even tunneling dynamics of Bose- 
Einstein condensates inside magneto-optical traps pj. In a condensed matter 
context, DNLS is usually derived from a coupled system consisting of an 
electron (or excitation), described quantum mechanically, that propagates 
in a discrete medium while interacting strongly with vibrational degrees of 
freedom, described by Einstein harmonic oscillators. In the limit where one 
assumes the oscillators completely enslaved to the electrons, one arrives at 
the DNLS equation: 

i (~^J = V ( C n+l + Cn-l) - X\Cn\ 2 C n (^ = 1), (!) 

where c n is the electronic probability amplitude at site n, V is the nearest- 
neighbor hopping parameter and \ is the nonlinearity parameter, propor- 
tional to the square of the electron- vibration coupling ||. The most remark- 
able feature of Eq.(|l|) is the existence of a critical nonlinearity parameter Xc, 
marking the onset of electronic selftrapping. The precise value of Xc depends 
upon the geometry and dimensionality of the lattice, but displays a rather 
universal behavior when properly scaled 0. 

Since most molecular degrees of freedom alway posses some degree of 
anharmonicity, it is interesting to examine how Eq.fll]) is modified when one 
considers weakly anharmonic and fully anharmonic oscillators. In this work, 
we investigate new DNLS-like equations obtained for three different choices 
of oscillator potentials: (i) weakly anharmonic: U(x) = (1/2) fcx 2 + (1/ 4)fc 3 x 4 , 
with k 3 <ti k (ii) Hard anharmonic: U(x) = k (cosh(x) — 1) and (iii) Soft 
anharmonic: U(x) = (k/3)[ u 2 + \u\ — log(l + \u\) ]. In particular, we focus 
on the case of a single vibrational impurity, and examine its bound state 
and dynamical selftrapping properties and compare them to the ones of the 
standard DNLS case. 
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1 Generalized DNLS Equations 



Let us consider an electron (or excitation) propagating in a onedimensional 
tight-binding lattice, and interacting with local vibrational degrees of free- 
dom modelled as identical, anharmonic Einstein oscillators at each site: 

i = v ( c n+i + Cn-i) + au n c n , (2) 

/ d 2 u n \ f dU\ 2 

where a is the electron-vibration coupling, m is the oscillator mass and U(u) 
is a general, anharmonic potential energy. In the m — > limit ("antiadia- 
batic" limit), the vibrations become completely enslaved to the electron and 
Eq.(|3|) becomes 

( dll \ 

-a \c n \ 2 . (4) 



du 



V(c n+1 + c n _i) + a u n (a\c n \ 2 ) c n , (5) 



After inverting Eq.(^), one obtains the oscillator equilibrium displacement 
in terms of the electronic probability, u n (t) = u n (a\c n \ 2 ) . Inserting this back 
into Eq.(||), leads to 

. / dc n s 
\dt / 

which constitutes formally a modified DNLS equation. Its specific form de- 
pends on the anharmonic potential chosen. In the special case of a harmonic 
potential U(x) = (l/2)kx 2 , © reduces to the well-known standard DNLS 
equation ([[]), with x = (ct 2 /k). 

2 Single Impurity 

In the case of a single vibrational impurity, located at the origin, Eq.(|^) 
simplifies to 

i (^fj = v (°n+i + Cn-i) + a WoO|c | 2 ) c 5 n fi- (6) 

We seek to determine the bound state(s) of Eq.(||). Taking c n = exp(-iEt) n , 
(H) becomes 

(E/V)(f) n = (pn+l + 4>n-l + (&/V) M (a|c | 2 ) 00 ^n,0- (7) 
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The simplest procedure now is to assume a solution of the form <p n = Ary n \ 
where A and 77 are real and \r)\ < 1. (Another, elegant method, is the use of 
lattice Green functions |7|, §]). At the impurity location, we have: (E/V) = 
2rj+ (a/V)u (aA 2 ) , while outside the impurity site, we have (E/V) = 77+1/77. 
After equating these two expressions, we obtain an equation for 77: 



rj 2 + 77 (a/V)u (aA 2 



1=0, 



with solution 



77 



" (w) MaA2) ± V 1 + (w) ul{aA2) 



(8) 



(9) 



Without loss of generality, we will take a > 0, (i.e., u (aA 2 ) < 0). In order 
to have \r)\ < 1, we have to take the negative root in (||). The unnormalized 
bound state profile then reads: 

|n| 



A 



a 



\2V 



u (aA z 



u 2 (aA 2 ) 



(10) 



The other equation needed comes from the normalization condition: 1 
X^oo I0n| 2 - This provides the following relation between 77 and A: 



V 



I - A 2 
1 + A 2 



(11) 



where rj is given by ([ 
equation for A: 

(1-A 2 \ 1/2 
[l + A 2 



After equating (|) and (|TTD , one obtains the following 



- (£) MaA') 



1/2 



(12) 



Once A (0 < A < 1) is determined, we can know uo(aA 2 ) and then <f> n , 
the bound state through ( |l0l) and the bound state energy E, from (E/V) = 
77 + (1/77), or 

f = -^l+(^) 2 ul(aA 2 ), (13) 

which always lies outside the band. Since we will consider oscillator potentials 
that increase monotonically with displacement, |Mo(av4 2 )| will always increase 
with increasing coupling causing, according to (0), an increasing detachment 
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of the bound state energy from the band, with increasing coupling. On the 
other hand, we can write the bound state as (fi n = A exp(— |n|/A), with 



A = - log 




-l 



(14) 



is the so-called "localization length". Then, as the coupling is increased, A 
decreases, giving rise to a more localized probability profile. 

3 Special Cases 

We consider now three, representative anharmonic potentials, characterized 
by U(u) monotonically increasing with \u\ and giving rise to restoring forces 
F(u) that are unbounded functions of u from above or below. This ensures 
that (|]) has always a real solution for the oscillator equilibrium displacement 
at the impurity site. In addition, we choose anharmonic potentials that 
are completely "soft" or "hard", that is, they always lie below or above 
the harmonic potential, for any displacement. With these conditions, we 
can predict the general features of the selftrapping behavior: Since U„ ^ < 
^harmonic < ^hard' ^ s means that the oscillator equilibrium displacement 
at the impurity site obeys: w^ aro - < ^harmonic < ^jjoft^ This implies that, 
the effective electron-vibration term in (P) is smaller (larger) for the hard 
(soft) case than for the harmonic one. In other words, the effective coupling 
is smaller (larger) for the hard (soft) case than for the usual, harmonic one. 

(1) Weakly anharmonic potential: U(u) = (l/2)ku 2 + (1/4)/c 3 m 4 . From 
(U) one obtains at the impurity site: kuo + k 3 u^ = —a A 2 , where A = <fto and 
we assume k 3 /k <C 1. Using a perturbative expansion, we have to first order 
in k 3 /k, 

V 



1 - 




(15) 



A simple graphical analysis of fll2|) , with uq given by (0), reveals the exis- 
tence of a solution for A, < A < 1, provided that a > \[2 (V — 1 = k). 
Figure 1 shows the electronic bound state probability profile \<j) n \ 2 for a small 
positive (negative) k 3 /k, corresponding to a weakly hard (weakly soft) un- 
derlying anharmonic oscillator. The harmonic (k 3 = 0) is also shown, for 
comparison. For a fixed coupling a, the bound state is wider (narrower) for 
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the hard (soft) case than for the purely harmonic case. This is to be ex- 
pected on general grounds, as explained before, since for a fixed coupling a, 
the effective coupling is smaller (larger) in the hard (soft) case than in the 
harmonic one. 

Using (|T5|) and the usual definition of the nonlinearity parameter \ = 
a 2 /k, the DNLS equation corresponding to this case is given by 

* = V (°n+i + c n-i) - x[ 1 - (h/k 2 )x\c \ 2 ] | c 1 2 Co S nt Q. (16) 

(2) Hard monotonically increasing potential: U(u) = /c(cosh(w) — 1). This 
is an example of a "hard" potential that is always greater than the harmonic 
one, approaching the harmonic case only in the limit of small displacement. 
From (S) one obtains: Mo — — sinh((a / k) A 2 ) . After inserting this into (|T2|), 
and performing some simplification, one arrives at the following nonlinear 
equation for A 2 : 

A 2 = ^VT^A* sinh- 1 ((a /k) A 4 ) (17) 



Simple analysis of ([L7|) reveals the existence of a single, nonzero solution, 
provided a > \[2 (V = I = k). The corresponding DNLS equation reads: 

i i^-^j = v (°n+i + c„-i) - a sinh -1 ((«/&) |C | 2 ) c 5 nfi . (18) 

(3) Soft monotonically increasing potential: U(u) = (k/3)[u 2 + \u\— log(l + 
\u\)}. This is an example of a "soft" potential that is always smaller than the 
harmonic one, approaching the harmonic case only in the u — > limit. From 



(D one obtains: u = (3/4) [ -(a/k) A 2 +1- ^1 + (2/3) (a/k) A 2 + (a/k) 2 A 4 }. 
After inserting this into (|1^) and after some simplification, the equation for 
A reads (V = 1 = k): 



A 2 = -aVT^A* 



aA 2 - 1 + J 1 + ^aA 2 + a 2 A 4 



(19) 



Now the situation is more complex, since a graphical examination of (|19|) 
reveals the existence of two critical coupling values = 1.36964 and = 
\/2, marking the boundary of three coupling regimes, with different number 
of bound states: For a < a)p , there is no bound state, at a = there is 
exactly a single bound state, with a finite amplitude. For a^' < a < af\ 
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there are two bound states. As a increases towards , one of the solutions 
approaches zero, while the other remains finite. Finally, at a = and 
beyond, there is a single bound state, whose localization length decreases 
with increasing coupling a. The DNLS equation for this case is 

dCn 



i l^J =V(c n+1 +c n . 1 )-a[(a/k)\c \ 2 -l + .Jl + (a/ky\Co\ i ]co6 nfi . (20) 

It is interesting to point out that the conditions for the existence of bound 
state(s) for our new DNLS equations ([18]) and (p0|) , differ significatively from 



the ones found in a previous investigation of bound states of a DNLS equation 
with a nonlinear term of the form x\^nf C n . There, the "soft" case < (3 < 2 
always displays a single bound state, while the "hard" case with (3 > 2, has 
a minimum nonlinearity threshold for \ below which no bound state exists 
and above which there are two bound statesjSJ. Now, as the reader can easily 
verify, this nonlinearity term \C n f can be thought of as originating from the 
(impurity) potential U = (k/j) |w| 7 , with 7 = 1 + (2/(3). Thus, (3 < 2 
corresponds to 7 > 2, that is, a hard case. Similarly, (3 > 2 corresponds to 
7 < 2, a soft case. In this light, the old results and the current ones look 
rather similar. 

Figure 2 shows the bound state energies for the hard, harmonic and soft 
cases, as a function of the coupling a. For the harmonic and hard cases the 
bound state energy curve starts at a = = \/2 and increases monotoni- 
cally in magnitude with increasing coupling. For the soft case, the eigenen- 
ergy curve starts earlier, at a = a)p = 1.36964. For coupling values between 
this and we have two eigenenergy curves, one of which moves away from 
the band upon increasing coupling, while the other approaches the band edge, 
reaching it exactly at a = as the inset in Fig. 2 shows. Figure 3 shows 
the bound state probability profile for the hard and soft cases, for a common 
coupling value a > Like the weak anharmonic case, the hard(soft) case 
has a larger (smaller) localization length than the usual harmonic case. 

We now compare the selftrapping transition for the three DNLS-like equa- 
tions ([TBI) , (IHf) and (f20"f). We place the electron initially completely on the 
impurity site (n = 0) and observe its time evolution. In particular, we focus 
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on the long-time average probability at the impurity site, 

(P)=lim (1/D [ T \C (t)\ 2 dt. (21) 

We use a fourth-order Runge-Kutta scheme whose precision is monitored 
through probability conservation J2 n \ c n(t)\ 2 = 1> an d a self-expanding lattice 
to avoid undesired end effects. Figure 4 shows (P) for the weakly anharmonic 
case, showing three curves corresponding to a weakly hard, harmonic and 
weakly soft oscillator. The standard selftrapping transition that occurs at 
a 2 /kV ~ 3.2 for k 3 = 0, is shifted to higher (lower) values for a small positive 
(negative) k%. This is consistent with the notion of an "effective" nonlinearity 
X e ff = x[ 1 — (A; 3 /A; 2 )x|co| 2 ] in (|i~6l). The selftrapping transition occurs when 
X e Q ~ 3.2. Finally, in Fig. 5 we show (P) for the fully anharmonic cases 
( |IBD and (pop, along with the standard (harmonic) case, for comparison. 
Like in the weak anharmonic case, the selftrapping transition is shifted to 
higher values for the hard case, while for the soft case it is shifted to a lower 
threshold. We also note that in this case, the selftrapping transition is more 
abrupt than in the other cases ||. 

4 Conclusion 

We have examined some new, DNLS-like equations that arise when consider- 
ing the propagation of an electron in a onedimensional lattice, while strongly 
interacting with local vibrational degrees of freedom modelled as anharmonic 
Einstein oscillators. We focused on the problem of a single vibrational im- 
purity, and computed the bound state(s) and the dynamical selftrapping 
properties. We found that, in general, the character of the anharmonicity 
affects the critical electron- vibrational coupling parameter needed to produce 
a bound state and to effect a selftrapping transition: For a hard (soft) anhar- 
monicity, the minimum coupling needed for a bound state to exist is larger 
(smaller) than for the usual, harmonic case. Similarly, the minimum nonlin- 
earity needed for selftrapping is larger (smaller) in the hard (soft) case than 
in the harmonic one. These results can be explained on the base of an "ef- 
fective" electron-vibration coupling that depends only on the magnitude of 
the oscillator equilibrium displacement at the impurity site, which in turns, 
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depends on the relative hardness of the anharmonic potential. 
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Captions List 



Fig.l : Weakly anharmonic vibrational impurity: Electronic bound state 
probability profile, for a small hard (k 3 > 0) and soft (k 3 < 0) anharmonicity. 
The harmonic case (k 3 = 0) is also shown, for comparison. (V — 1, k — 1, a — 
2 and k 3 /k = 0.2). 

Fig. 2 : Fully anharmonic vibrational impurity: Electronic bound state en- 
ergy as a function of the electron-vibration coupling, for the fully hard and 
soft cases considered in section 3 (V — l,k — 1). The usual harmonic case 
is also shown, for comparison. 

Fig. 3 : Fully anharmonic vibrational impurity: Electronic bound state 
probability profile for the fully hard and soft cases considered in section 3 
(V = l,k = 1 and a = 2). 

Fig. 4 : Weakly anharmonic vibrational impurity: Long-term average prob- 
ability for finding the electron at the vibrational impurity site, as a function 
of the nonlinearity parameter x — ® 2 /k. 

Fig. 5: Same as in Fig. 4, but for the fully anharmonic case. 
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